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Abstract 

The QCD improved parton model is a very successful concept to treat processes in hadronic 
interactions, whenever large partonic transverse momenta p± are involved. However, cross sec- 
tions diverge in the limit p± — > 0, and the usual treatment is the definition of a lower cutoff 
P_i_min, such that processes with a smaller p± - so-called soft processes - are simply ignored, which 
is certainly not correct for example at RHIC energies. A more consistent procedure amounts 
to introduce a technical parameter Qq, referred to as soft virtuality scale, which is nothing but 
an artificial borderline between soft and hard physics. We will discuss such a formalism, which 
coincides with the improved parton model for high p± processes and with the phenomenological 
treatment of soft scattering, when only small virtualities are involved. The most important 
aspect of our approach is that it allows to obtain a smooth transition between soft and hard 
scattering, and therefore no artificial dependence on a cutoff parameter should appear. 



1 Introduction 

The standard parton model approach to hadron-hadron scattering amounts to presenting the partons 
of projectile and target by momentum distribution functions, /j and fj, and calculating inclusive cross 
sections for the production of parton jets with the squared transverse momentum p\ larger than some 
cutoff pimin as 

cr- m c(s)=J2 dp 2 ± dx + dx'fi(x + ,p\)f j (x~,p\)da ij /dp 2 L (x + x's) (1) 

ij J P 2 ±>p 2 ±mixl 

where d&ij/dp\ is the elementary parton-parton cross section and i,j represent parton flavors. Many 
Monte Carlo applications are based on the above formula: ISA JET gj, PYTHIA @,g|,HERWIG @-|], 
HIJING - for more details see the review ||. 
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It is the cutoff p\ min which prevents the above integral to diverge, and there are several justifica- 
tions for using such a cutoff. If one is only interested in hard processes, for example to study large 
Pi_ jets, such a cutoff can be used without problems, since soft reactions are simply not considered. 

The situation is much more complicated when it comes to treating an average event, since here 
the soft component plays an important role. In nuclear scattering, one may argue with screening 
effects: in heavy nuclei at very high energies, the parton density will be reduced due to screening, and 
this may be effectively done via a cutoff. Real life is more complicated: we have finite nuclei at finite 
energy, so soft physics still plays a role and has to be treated properly. The same is true for screening 
corrections: in particular for Monte Carlo applications one cannot treat them just introducing a 
cutoff or modifying structure functions, this is only valid for calculating inclusive spectra based on 
the assumption that factorization holds. 

There is already some indication that the above-mentioned concept comes to its limits: taking 
parton distribution functions compatible with latest HERA measurements, it seems to be impossible 
to get the correct energy dependence of the proton-proton cross section, which can only be cured by 
an unreasonable assumption like an energy dependent cutoff (see, for example, 0). 

We think that one should well separate the aspect soft/hard physics and screening. One may 
introduce a soft scale Qq, which is meant to be the borderline between soft and hard physics, such 
that above this scale one may use perturbative methods, whereas below one has to work with param- 
eterizations. We do not think that the physics changes abruptly when crossing this artificial border, 
it is more a technical problem that we know how to do calculations only above the scale Qq. It is 
not at all a serious problem that one has to rely on a phenomenological treatment at low virtualities, 
since there exists a wealth of information about this topic, in particular since such concepts as the 
Pomeron can be investigated at the HERA collider. Being just a borderline between soft and hard 
physics, there is absolutely no reason that the parameter Ql should be energy dependent, or vise 
versa - introducing an energy dependent cutoff necessarily implies changing physics content of the 
"soft" part of the interaction as with the increasing energy harder and harder partons are treated as 
being "soft". 

Naturally, at high energies the role of screening corrections due to so-called enhanced diagrams 
becomes extremely important. Our approach allows to treat corresponding contributions explicitely, 
as coming from Pomeron-Pomeron interactions in the "soft" nonperturbative region. This allows to 
achieve a microscopic description of the interaction process, to resolve the seeming contradiction 
between the energy dependence of hadronic interaction cross sections and the realistic structure 
functions measured by HERA, and to get rid of the artificial dependence of the results on the 
technical parameter Ql ||T0| . 

2 Qualitative Discussion of Nucleon-Nucleon Scattering 
2.1 The Structure of the Nucleon 

Nucleons are composite objects, and high energy nucleon-nucleon scattering requires therefore some 
discussion on the intrinsic structure of the nucleon. Let us consider deep inelastic scattering, where a 
virtual photon probes the nucleon structure, see fig. |l], where we show a cut diagram, corresponding 
to a squared amplitude. As shown in the figure, the photon couples to a quark from the nucleon, 
allowing this way the measurement of the quark momentum distributions (nucleon structure func- 
tions). For photons with large virtualities, one may apply the methods of perturbative QCD in 
order to understand the structure of the "blob" in the above figure: we know that the final quark 
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Figure 1: A virtual photon probing the structure of the nucleon. 

is emitted from a parton, which itself is emitted from another parton and so on. Employing the 
leading logarithmic approximation one obtains the following picture (see fig. ||): The nucleon emits 




Figure 2: A highly virtual photon couples to the "end" of a parton ladder, emitted from the nucleon. 

a parton with a virtuality of the order of some soft scale Ql and with some light cone momentum 
fraction x, which emits itself parton with higher virtuality but smaller x and so on, such that the 
last parton couples to the photon. It is a question of taste to consider the parton ladder of being a 
part of the internal structure of the nucleon and the photon to interact with a quark of the nucleon, 
or to consider only the first parton at scale Ql (the "soft" parton) being a constituent of the nucleon, 
whereas the parton ladder is considered to be a part of the interaction. Concerning the soft partons, 
one distinguishes between sea quarks or gluons and valence quarks. Whereas the latter ones have 
typical light cone momentum fraction distributions of the form x~ 5 , the former ones are distributed 
as x~ a °, with «o being at least unity. This means that sea quarks or gluons have typically very 
small x-values, contrary to the valence quarks. What is the mass of the "blob" in fig. [|, left behind 
the emitted quark? Having emitted a parton with momentum k and virtuality k 2 = —Q\ from the 
nucleon, considered forward moving with momentum p, we obtain in the high energy limit for the 
squared mass of the remainder (nucleon minus parton) 

s = -p+k~ = p + Q 2 Jk + = Ql/x, 

where p ± = p ±p 3 , = k ±k 3 are the light cone momenta of the nucleon and of the emitted parton 
correspondingly and x = k + /p + is the light cone momentum fraction of the parton. This shows that 
in case of small x (sea quarks or gluons), the remainder has a large mass; whereas for valence quarks, 
we obtain masses of the order of Qq. This means sea quarks or gluons are not elementary nucleon 
constituents, but they are rather emitted from some object with large mass and low virtuality - 
Such objects are usually referred to as soft Pomerons, and very powerful phenomenological techniques 
exist to deal with them. 

So we distinguish between two contributions, depending on whether the first (soft) parton is a 
sea quark or gluon or a valence quark. In case of valence quarks, we consider the quark to be a 
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direct constituent of the nucleon, and the corresponding photon-nucleon scattering diagram is shown 
in fig. |3|. So one may assume some initial distribution F va \(x) for the valence quarks, and treat 
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Figure 3: Photon-nucleon scattering, with the first (soft) parton being a valence quark. 

then the scattering of this valence quark with the virtual photon. In case of sea quarks or gluons, 
there has to be a soft Pomeron between the nucleon and the first parton, where we do not know the 
precise microscopic structure, but we know how to parameterize the corresponding amplitudes. So 
our graphical representation will be simply a "blob", as shown in fig. £|. This blob stands for a soft 




Figure 4: Photon-nucleon scattering, with the first (soft) parton being a sea-quark or gluon. 

Pomeron, and we will use later the corresponding parameterizations. We do not know how the soft 
Pomeron couples to the nucleon, but we may imagine "some constituent" carrying some fraction x 
of the nucleons momentum according to some distribution F(x), which emits a soft Pomeron, which 
by itself emits the first perturbative (k 2 ~ — Ql) parton. In this sense, the external lower line in fig. 
f| may be considered to be such a constituent. 

To summarize: we consider in any case the scattering of a virtual photon with a nucleon con- 
stituent carrying some fraction x of the nucleons momentum. In case of valence quarks, this con- 
stituent is simply a valence quark with distribution F va i(x), in case of sea quarks or gluons, the 
constituent is not known, and we use some function F(x), to be specified later. In case of valence 
quarks, the interaction process is treated entirely in the framework of perturbative QCD, using the 
leading logarithmic approximation; in case of sea quarks or gluons, one has in addition to consider a 
soft Pomeron. 



2.2 Nucleon-Nucleon Scattering 

Virtual photon-nucleon scattering may be seen as an interaction between the virtual photon and a 
nucleon constituent, where the interaction is realized via the exchange of some composite object - as 
discussed in the preceding section. This exchange object is either a parton ladder (in case of valence 
quarks involved) or a parton ladder plus a soft Pomeron (in case of sea quarks or gluons involved). We 
easily generalize to nucleon-nucleon scattering. Let us assume that there are hard partons involved 
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- with virtualities bigger than the soft scale Qq. Employing the leading logarithmic approximation 
of perturbative QCD, we expect some hard scattering of two partons with high virtualities, these 
two parton being emitted from partons with smaller virtualities and so on, till one reaches on both 
ends the soft scale Qq. So now we have two initial (low virtuality) partons on both sides, and for 
either of them we can literally repeat the discussion of the last chapter: each parton may be of 
sea or valence type, and correspondingly we have four contributions. If both partons are of valence 
type, we have an interaction between two valence quarks via the exchange of a parton ladder, as 
shown in fig. |5] (valence- valence contribution). We do not show in the figure the spectator quarks. 
The fact that the ladder rungs getting narrower towards the middle of ladder indicates symbolically 
virtuality ordering: a bigger ladder rung represents smaller virtuality. Each of the two ordered ladders 
corresponds exactly to the parton ladder employed in deep inelastic scattering. The only new aspect 




Figure 5: The valence-valence contribution. 

here is the fact that the two ordered ladders are "glued" together by means of an elementary Born 
scattering diagram, which represents the process with the highest virtuality involved. Since only 
hard partons contribute, we refer to the valence-valence contributions also as the "hard" one. 

We obtain a second contribution, when we consider the scattering of sea quarks or gluons, as 
shown in fig. |. As discussed in the preceding section, there is in general a soft, large mass object 




Figure 6: The sea-sea contribution. 

between the sea quark (gluon) and the nucleon, which we identify with the soft Pomeron. So the 
sea quarks or gluons are internal lines in the above diagram. The external lines in the figure refer to 
"nucleon constituents", the precise nature of which we do not need to specify at this point. We simply 
assume that they carry a fraction x of the nucleons momentum, according to some distribution F(x). 
We have of course also the mixed cases, where a valence quark scatters off a sea quark or gluon, as 
shown in fig. |7|. The external legs here are a valence quark on one side and a nucleon constituent on 
the other side, since the sea quark or gluon is as usual an internal line. 
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Figure 7: Two "mixed" contributions. 

Let us come back to the sea-sea contribution. Here we consider hard perturbative partons in 
the ladder and "soft partons" in the shaded area before it, where the latter ones are not treated 
explicitly but rather parameterized as soft Pomeron contributions. One may imagine the case where 
the hardest parton in the process is already a soft one, with the corresponding virtuality being smaller 
than Qq cutoff, so the hard piece in the middle is reduced to zero. Since this is not possible for our 
sea-sea contribution, where one requires always at least one hard parton in the diagram, one has 
to add explicitly a purely soft contribution, which is simply parameterized as a soft Pomeron and 
graphically represented as a "blob", see fig. |[ We use of course the same parameterization as for 




soft 

Pomeron 



Figure 8: The soft contribution. 

the soft Pomerons appearing in the sea-sea and the mixed contributions, as we are going to discuss 
later. The external legs for the soft contribution are nucleon constituents, as for the 

A complete elementary interaction in hadron-hadron scattering is therefore the sum of all the 
above-mentioned contributions: a soft one, a hard one, and three "semi-hard" contributions (sea-sea 
and mixed). We have a smooth transition from hard to soft physics: lowering the energy will reduce 
the relative weight of the semi-hard contributions till they finally die out and we are left with a 
purely soft contribution. In our approach, we do not consider soft and hard physics as completely 
different: the soft domain is just the continuation of perturbative domain into a region, where we do 
not have the technical abilities to perform rigorous calculations. But physics changes smoothly. 

In the following section, we are going to treat the different contributions in detail. 

3 Detailed Treatment of Nucleon-Nucleon Scattering 
3.1 Soft Scattering 

Let us first consider the purely soft contribution of the fig. || where all virtual partons appearing 
in the internal structure of the diagram have restricted virtualities Q 2 < Qq, where Ql ~ 1 GeV 2 
is a reasonable cutoff for perturbative QCD being applicable. Such soft non-perturbative dynamics 
is known to dominate hadron-hadron interactions at not too high energies. Lacking methods to 
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calculate this contribution from first principles, it is simply parameterized and graphically represented 
as a 'blob', see fig. |[ It is traditionally assumed to correspond to multi-peripheral production of 
partons (and final hadrons) jl2| and is described by the phenomenological soft Pomeron exchange 
contribution. The scattering amplitude is given as (TT 



T so ft(M) = 87rs r/(t) 7^ art f — J exp(A2 t (s/s )t) , (2) 

with 

^lol( Z ) = ^Jart + "soft ln^, (3) 

where s = (p + p') 2 , p,p' are the 4-momenta of the initial constituent partons in the process before 
and after the scattering (considered here as being nearly real, i.e. p 2 = p' 2 ~ ). The parameters 
av>ft(0), a'soft are the intercept and the slope of the Pomeron trajectory, 7 part and -Rp art are the vertex 
value and the slope for the Pomeron-parton coupling, and sq ~ 1 GeV 2 is the characteristic hadronic 
mass scale. The so-called signature factor r\ is given as 

v (t) = i - cot — ^ ~ i. (4) 

Cutting the diagram corresponds to the summation over multi-peripheral intermediate hadronic 
states, connected via unitarity to the imaginary part of the amplitude (||), 

^disc^ T soft (s, t) = ^[T so f t (s+iO,t) -T soft (s-iO,t)} (5) 
= 2ImT soft (S,t) (6) 



n, spins,. 



\ dT nTp,p'^X n Tlp>^X n i ( 7 ) 



where T PiP i^x n is the amplitude for the transition of the initial partons p, p' into the n-particle state 
X n , dr n is the invariant phase space volume for the n-particle state X n and the summation is done 
over the number of particles n and over their spins and species, the averaging over initial parton 
colors and spins is assumed; disc^ T so ft(s, t) denotes the discontinuity of the amplitude T so f t (s,t) on 
the right-hand cut in the variable s. 

The corresponding profile function for parton-parton interaction, defined as twice the imaginary 
part of the Fourier transform T of T, divided by the initial parton flux 2s, 

D soit (s,b) = — 2Imf soft (s,6), (8) 



is given as 



1 ( 

D soft (s,b) = jd 2 q±exp(-iqj?) 2ImT soft (s, -q 2 ^j (9) 

(10) 



27 2 



part 




For t = one obtains via the optical theorem the contribution cr so f t of the soft Pomeron exchange to 
the total parton interaction cross section, 

^soft(s) = ^2ImT soft (s,0) (11) 
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d 2 bD soit (s,b) (12) 

(A \ Oisoft(O)— 1 
v ■ (13) 

The external legs of the diagram of fig. || are "partonic constituents", as discussed in the preceding 
section. 

3.2 Hard Scattering 

Let us now consider the hard (or valence-valence) contribution, see fig|| All the processes in the 
interaction diagram are perturbative, i.e. all internal intermediate partons are characterized by 
large virtualities Q 2 > Qq. In that case, the corresponding hard parton-parton scattering amplitude 
Ih ard (s, t) (J, k denote the types (flavors) of the initial quarks) can be calculated using the perturbative 
QCD techniques fll4], fi7j||, and the intermediate states contributing to the absorptive part of the 



amplitude can be defined in the parton basis. In the leading logarithmic approximation of QCD, 
summing up terms where each (small) running QCD coupling constant a s {Q 2 ) appears together with 
a large logarithm ln(Q 2 /Aq CD ) (with Xqcd being the infrared QCD scale), and making use of the 
factorization hypothesis, one obtains the contribution of the corresponding cut diagram for t = q 2 = 
as the cut parton ladder cross section <J^ Td (s, Ql) [], which will correspond to the cut diagram, where 
all horizontal rungs are the final (on-shell) partons and the virtualities of the virtual t-channel partons 
increase from the ends of the ladder towards the largest momentum transfer parton-parton process 
(indicated symbolically by the 'blob' in the middle of the ladder): 



<Us,Ql) = ±21mT&l$,t = 0) (14) 
^cd04, Ql M 2 F ) i?S CD fe, Ql M 2 F )6(M 2 F - Ql) . (15) 



ml 

X ' " ' 



Here da^ vn / dp\ is the differential 2 — > 2 parton scattering cross section, -p\ is the parton transverse 
momentum in the hard process, m, I and x% are correspondingly the types and the shares of the 
light cone momenta of the partons participating in the hard process, and Mp is the factorization 
scale for the process (we use Mp = Pj_/4). The 'evolution function' Eq^, d (QI, Mp, z) represents the 
evolution of a parton cascade from the scale Ql to Mp, i.e. it gives the number density of partons 
of type m with the momentum share z at the virtuality scale Mp, resulted from the evolution of the 
initial parton j, taken at the virtuality scale Qq. The evolution function satisfies the usual DGLAP 
equation [16|] with the initial condition Eq^ d (Qq, Ql, z) = 5 J m 5(1 — z). The factor K ~ 1.5 takes 
effectively into account higher order QCD corrections. 

In the following we shall need to know the contribution of the uncut parton ladder T^ rA (s,t) 
with some momentum transfer q along the ladder (with t = q 2 ). The behavior of the corresponding 
amplitudes was studied in |17|] in the leading logarithmic(l/x) approximation of QCD. The precise 
form of the corresponding amplitude is not important for our application; we just use some of the 
results of |17|], namely that one can neglect the real part of this amplitude and that it is nearly 
independent on t, i.e. that the slope of the hard interaction i?^ ard is negligibly small, i.e. compared 

1 Strictly speaking, one obtains the ladder representation for the process only using axial gauge. 
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to the soft Pomeron slope one has -Rh ard — 0- So we parameterize T^ rd (s,t) in the region of small t 



as 



18 



T hard(s, t) = is a^ rd (s, Ql) exp (r*^ t) (16) 

The corresponding profile function is obtained by calculating the Fourier transform Th ar d of Th ar d 
and dividing by the initial parton flux 2s, 



which gives 



Dil d (S,b) = ±21mf£ d (s,b), (17) 



DliA^b) = ^Jd 2 q ± e W (-iq ± b)2ImT^ ld (s,-q 2 ± ) (18) 



°hard 



^^lM~^lJ (19) 



Since we also talk about "valence-valence" contribution, we use sometimes -D va i-vai instead of 

-Dhard^ 

£>S-vai {s,b) = DH d (s,b) , (20) 
so these are two names for one and the same object. 

3.3 Semi-Hard Scattering 

Let us start discussing semi-hard sea-sea contribution, represented by a parton ladder with "soft 
ends", see fig. ||. As in case of soft scattering, the external legs are nucleon constituents, connected to 
soft Pomerons. The outer partons of the ladder are on both sides sea quarks or gluons (therefore the 
index "sea-sea"). The central part is exactly the hard scattering considered in the preceding section. 
With a sum over all the hard scattering processes in the center, we get the mathematical expression 
for the corresponding amplitude 

^T sea _ sea (s,t)=£ f iT^ ld (z + z~s,t) JmT^J^t) ImT s fc of /^,t) , (21) 

j k 

with z being the momentum fraction of the external legs of the parton ladder relative to the 
momenta of the nucleon constituents. The indices j and k refer to the flavor of these external ladder 
partons. The amplitudes T s J oft are the soft Pomeron amplitudes discussed earlier, but with modified 
couplings, since the Pomerons are now connected to the ladder on one side. The arguments Sq/z ± 
are the squared masses of the two soft Pomerons, z + z~s is the squared mass of the hard piece. 

Performing as usual the Fourier transform to the impact parameter representation and dividing 
by 2s, we obtain the profile function 

Aoa-sea (s, 6) = 7^ 21m f sea _ sea (s, b) , (22) 

which may be written as 

Aca-sea (a, 6) = £ Cdz+ dz~E{ oit (z + ) E^ o{t (z~) at d {z + z~s, Q 2 Q ) (23) 

jk J0 

1 / b 2 \ 

x , :™ ; 7. . ^ ex P 



4^A ( s 2 ) ft (l/(z+^)) V 4^(1/0*+*-)). 



9 



with the soft Pomeron slope A soft and the cross section a 3 hwd being defined earlier. The functions 
i?soft (z ± ) representing the "soft ends" are defined as 

Ho^ ± )=^TU^,t = 0). (24) 



With the phenomenological description of soft Pomeron exchange ||13|| , we parameterize the amplitude 



soft aS 



T s s oft (s, t) = 87ts oV (t) 7part79 (-T* <°> exp (A« (s/s )t) (1 - (25) 

So s 



^oftW=^a rt + «softln^ (26) 

E 9 soit (z) = 8ns 0l ^a 9 ^ Qsoft(0) (1 - zf' . (27) 



32*(*)=7« / d£P«(OE a soit (~) , (28) 



where 

So we get £ s g oft : 
With (|27|), we obtain 

with 

7<? 9 7 S = Wgpfo %, Ig = (1 - ^split) 7 S , (29) 
where tf S piit parameter determined the relative weight of sea quark content of the soft Pomeron (the 
probability for g qq splitting in the soft Pomeron) . By the momentum conservation constraint 

1 = f 1 dz^z Host (*) = 87r So7p a rt 7 9 / dz z 1 -*^ (1 - zf« , (30) 

j 

7 9 is fixed as 

~ = 1 T(3-a soi (0)+f3 a ) (31) 

9 87r S o 7 p art r(2-a soft (0)) r(l + /3 g )' 1 ' 

We neglected the small hard scattering slope -Rh ard compared to the Pomeron slope A so ft- We call 
Egoft also the " soft evolution", to indicate that we consider this as simply a continuation of the QCD 
evolution, however, in a region where perturbative techniques do not apply any more. It's easy to 
see that -Eg 0ft (z) has the meaning of the momentum distribution of parton j in the soft Pomeron at 
virtuality Qfo. 

Consistency requires to also consider the mixed semi-hard contributions with a valence quark on 
one side and a sea quark or gluon on the other one, see fig. |7|. We have 

<i-sea(M) = jf^EM^^.t) (z-8,t) (32) 

J Z j, \ Z J 

and 

BL-M-X/y^) < d (,-i,g S ) -^ r (--J^ (33) 

where j is the flavor of the valence quark at the upper end of the ladder and k is the type of the 
parton on the lower ladder end. Again, we neglected the hard scattering slope -R^ard compared to 
the soft Pomeron slope. A contribution -Df ea _ va i (s,b), corresponding to a valence quark participant 
from the target hadron, is given by the same expression, 

^ea-val & 6) = ^al-sea & &) , (34) 

since eq. (|33|) stays unchanged under replacement z" — > z + and only depends on the total cm. 
energy squared s for the parton-parton system. 



10 



3.4 Nucleon-Nucleon Scattering 

Let us define for any of the elementary interaction types K (soft, hard, val, and val-sea) 

a "dressed" partonic profile function G via 

G K {x + ,x~,s,b) = — 2lmf' K (x + ,x~,s, b), (35) 

2x + x~s 

where T'k is the Fourier transform of T'k, with 

T' K (x + ,x-,s,q 2 )=T K (x + x-s,q 2 ) F K+ {x + )F K -{x-) exp(2R 2 N q 2 ) (36) 

representing the contributions of "elementary interactions plus external legs". Here we assumed a 
simple factorized form for the vertex for nucleon coupling to n constituent "soft" partons participating 
in elementary scattering processes (external partons for processes of fig. [|, H, 0, H) : 



N ( n\p, h,..., k n , q x , . . . , q n ) = J| 



n 

2 

i 

i=lL \ k=l 



I 

F Kt (x±)exp(R 2 N J2<ll 




F Temn \l-^xf\ , (37) 



where x ± = ^/p^ are the light cone momentum fractions of the i-th constituent parton (± corre- 
spond to the projectile/target case), g« is the 4-momentum transfer in the i-th scattering process, Kf 
is the type of i-th constituent parton (Kf ="sea" for the "soft" or "sea-sea" type process; Kf ="val" 
for the "hard" scattering; Kf ="val", K~ ="sea" for the "val-sea" case etc.), and R 2 N is known as 
the nucleon Regge slope. The factors in the square brackets are then associated with individual 
scattering contributions and included in the definition of T' k ,Gk- 

For the functions F sea , F remn we use a simple Regge-inspired ansatz: 

F sea (x) = lN x~ a ^\ (38) 

F Temn (x) = x a — . (39) 

where the different parameters a are finally determined by comparing with experimental data. The 
function F va i is constructed such that it reproduces an input parametrization (GRV94 |19|) for valence 
quark momentum distributions in the nucleon q v (x, Qq) at the initial scale Qq . 
Based on the above definitions, we may write 

G K (x + ,x~,s,b) = D' K (x + aTs,&) F K +(x + ) F K -(x~), (40) 

with the .D'-functions expressed in terms of the bare parton profile functions as 



8R% 



(41) 



which means that D'k has the same functional form as Dk, with X^ n \^) being replaced by A^(£) + 
2R 2 N . 

So we have the "bare" partonic profile functions Dk and the "dressed" ones Gk, where the former 
ones describe pure parton-parton |20| scattering (or more generally scattering between partonic 
constituents), whereas the latter ones describe partonic scattering, however, taking into account the 
longitudinal and transverse momentum distribution of the partons in the nucleons. These dressed 
profile functions are the elementary quantities, based on which one may construct a proper multiple 
scattering theory, for nucleon-nucleon as well as nucleus-nucleus scattering. 
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Before we discuss a more correct treatment, let us first show how we recover the conventional 
Gribov-Regge approach of multiple scattering. Here one considers so-called eikonals, 

Xk(s, b) = J dx + dx~G K (x + , x~ , s, b) F remn (l F remn (l -ar~) , (42) 

and using the eikonal approximation, one may obtain the inelastic cross section as 



0"inel(S 



J d 2 b jl — exp 



Y,XK(s,b) 

K 



(43) 



To obtain such a simple formula, one has to ignore the fact that the energy has to be shared between 
the different elementary interactions, which is certainly not correct. 
Doing energy conservation properly, one obtains ||20| 

rf 2 6^— / \{dx+dx- [] G(x+,x~,s,b) E 77 / Y[dx$dx x 
^1 m\ Jo ^ » ^ ^ l\ Jo 

x II -G(xt,%x, s, b) F rcmn ^ proj -X: xt^j F rcmn (x^-J^ x^j , (44) 



with x^,x^ being the light cone momentum fractions of constituent partons for the elementary 
scattering contributions and 

1"E< (45) 



x 



proj/targ 



being the momentum fraction of the projectile/target remnant, and with 

G G so ffc -\- G^val — val G sc£l — sea -|- (jr va j_ sea -|- G sea _ va j (46) 

being the sum of all the contributions, soft ones, hard ones, and semi-hard ones. The terms +G 
in formula (]44D correspond to cut elementary diagrams, the terms — G to the uncut ones, the sums 
represent summations over all possible cuts. Eq. fl44D can be written as 



<x ine i(s) = / d 2 6{l-$(s,6)}, (47) 

with 



®(s,b)=itjj I J[dxtdx x H -G(x$,x x ,s,b)F iemn (l x$ \ F rcmn 

1=0 L J ° A=l A=l V A / V A 



A > 



(48) 



which can be evaluated numerically. This energy sharing formalism can be easily generalized to 
nucleus- nucleus scattering [BO], not only for calculating cross sections, but also particle production. 



4 Some Results and Discussion 

4.1 The Behavoir of the "Dressed Profile Function" G 

As discussed in detail in the preceding chapter, one may express the inelastic cross section (and 
many other quantities) in terms of the "dressed profile function" G, which is a sum of soft, hard, 
and semi-hard contributions: G = J2Gk- This function can be interpreted as the number of ele- 
mentary interactions ("Pomerons") with light-cone momentum fractions x + and x~ at a given impact 
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Figure 9: The functions Gk in dependence oni = x + x~sX 6 = for the different contributions K for a cms energy 
of 200 GeV. We show soft (dotted), semi-hard (dashed), and valence contributions (dashed-dotted), where semi-hard 
is meant to be the sum of sea-sea, sea-val, and val-sea. 
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Figure 10: The functions Gk in dependence oni = x + x~aX b — for the different contributions K for a cms energy 
of 1800 GeV. We show soft (dotted), semi-hard (dashed), and valence contributions (dashed-dotted), where semi-hard 
is meant to be the sum of sea-sea, sea-val, and val-sea. 



parameter b. We first investigate the functional form of these functions Gk- We plot in figs. || and 
|10| the functions Gk in dependence on x = x + x~at b = for a cms energy of 200 and 1800 GeV 
for the different contributions K (we consider the sum of sea-sea, sea-val, and val-sea, referred to 
as "semi-hard"). The dominant contributions are in any case the soft and the semi-hard ones. Both 
show roughly a power-law behavior, with a much steeper increase of the semi-hard component. At 
200 GeV, the soft component is by far dominant, whereas for 1800 GeV the semi-hard component is 
taking over, however, with the soft one still being bigger at small x. Although the semi-hard part 
has a cutoff at some small value of x, the complete contribution shows always a smooth and regular 
behavior in the limit of x going to zero, due to the soft component. 

The smooth behavior of the G function has also the very nice side effect that we can approximate 
it (to a very good precision) by a simple analytical function, which allows an analytical calculation of 
many interesting quantities, in particular in nucleus-nucleus collisions, where numerical calculations 
are quite costly, or even impossible. 

Having calculated G, we can now proceed to calculate the inelastic cross section, using eqs. 
(|47|J48D . We obtain the result shown in fig. || as a full curve. Accounting for lowest order screening 
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Figure 11: The inelastic cross section for proton-proton scattering as a function of the energy y/s. 



corrections due to so-called enhanced Pomeron diagrams we obtained somewhat reduced values of 



0"inei(s), shown in fig. |TT] as the dashed curve. 



4.2 Comparison with Data: 

An important advantage of our approach is that it allows not only to calculate different interaction 
cross sections but also to develop a fully self-consistent Monte Carlo procedure to simulate hadron- 
hadron (nucleus-nucleus) interactions, including direct modeling of perturbative parton evolution in 
hard and semihard elementary processes. Thus, we were able to simulate many processes and obtain 
valuable results which agree quite well with experimental data. Here we present our results for 
proton-proton interactions in the energy range between roughly 10 and 2000 GeV, which represents 
the range of validity of our approach. The energy dependence of the total cross section has been 
compared with the data in figure [12]. The average multiplicities of different hadron species as a 



function of the energy y/s are shown in figure [13 
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Figure 12: The total cross section as a function of the energy y/s (left figure): the full line is the simulation, the 
points represent data j2l]]. Pomeron numbers as a function of energy y/s (right): soft (dashed) and semi-hard (solid 
line) Pomerons. 
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Figure 13: The average multiplicities of different hadron species, as a function of the energy y/s. From top to bottom: 
all charged particles, 7r + , 7t~, K + ,K~ jj. The full lines are simulations, the points represent dada from (23]. 
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Figure 14: Rapidity distributions of pions at 100 GeV. The full lines are simulations, the points represent data |2J 



We also present the spectrums of pion here. Figure [14] shows the rapidity distributions of pions at 
100 GeV and figure [15| shows the psedo-rapidity distributions of pions (7r + ,7r~)and charged particles 
(all charged and negatively charged) at 200 GeV. Figure [16] shows the psedo-rapidity distributions 
of positively and negatively charged particles at 53 GeV (cms). Figure [17| shows the psedo-rapidity 
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Figure 15: Psedo-rapidity distributions of pions (ir + ,ir )and charged particles (all charged and negatively charged) 
at 200 GeV. The full lines are simulations, the points represent data pa,E4]. 
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Figure 16: Psedo-rapidity distributions of positively and negatively charged particles at 53 GeV (cms). The full lines 
are simulations, the points represent data pa]. 



distributions of charged particles at 200 and 1800 GeV (cms). 

The spectrums of protons are presented, too. Figure |l8| shows the longitudinal momentum fraction 
distribution of protons at 100 GeV, integrated over p t (left) and for p t = 0.75 GeV/c (right). Figure 
|l~9l shows the longitudinal momentum fraction distribution of protons at 200 GeV, integrated over p t 
(left) and at 175 GeV for p t = 0.75 GeV/c (right). 
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Figure 17: Psedo-rapidity distributions of charged particles at 200 and 1800 GeV (cms). The full lines are simulations, 
the points represent data ]2^,[2^|. 
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Figure 18: Longitudinal momentum fraction distribution of protons at 100 GeV, integrated over p t (left) and for 
pt = 0.75 GeV/c (right). The full lines are simulations, the points represent data p8| . 
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Figure 19: Longitudinal momentum fraction distribution of protons at 200 GeV, integrated over p t (left) and at 175 
GeV for pt = 0.75 GeV/c (right). The full lines are simulations, the points represent data |23|,|2"s|]. 
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4.3 Some Discussion 

This paper has shown a theoretical discussion of a new concept of treating soft and hard scattering 
in a consistent fashion. As discussed above, at energies presently (or in the near future) accessible, 
both soft and hard components play an important role, so one cannot simply ignore any of the two. 
This was our main motivation to develop the theoretical framework discussed in this paper. The full 
self-consistent scheme for the description of high energy hadron-hadron (nucleus-nucleus) interactions 
has to account also for all significant screening corrections, or more generally interactions between 
elementary diagrams. This part of the work will be discussed elsewhere UToH - 
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